Geodesies on deformed spheres asymptotically described using the Funk transform 
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We consider geodesies on the surfaces obtained by weak deformations of the standard 2D-sphere. 
The dynamics of a particle on the surface can be asymptotically described by the averaged evolution 
of the particle's angular momentum. It is shown that the system describing this evolution has a 
Hamiltonian, which is obtained by applying the Funk transform to the function defining the deviation 
of the surface from the standard sphere. This system has the 2D-sphere as its phase space, so it is 
integrable and its trajectories admit of topological description in terms of its phase portrait on the 
sphere. 
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I. INTRODUCTION 



i-C . The problem of finding geodesies on a surface is a classical subject of differential geometry and analytical mechanics. 

Generally the exact analytical treatment of the dynamics equations gets quite involved even for comparatively simple 
surfaces such as ellipsoid (the well-known result of Jacoby, On the other hand, if we are interested in the 



qualitative topological structure of the system's trajectories and their large-scale properties, exact solutions are not 
necessary. For such purposes we can apply asymptotic methods, which provide a general picture of the system's 
dynamics. 



In the present work we consider surfaces that are obtained by a small deformation of the standard 2D-sphere. We 
■ perform an asymptotic Hamiltonian reduction of the dynamics of a particle on such a surface to a Hamiltonian system 
Q\ ' with one degree of freedom. The latter system is always integrable and admits of visual description by means of phase 
portraits. The idea underlying the reduction was introduced in Q, where a special class of deformed spheres was 
analyzed in detail. The key observation is that for small periods of time the trajectory of the particle is close to the 
exact solution of the unperturbed problem, i.e. a geodesic on the standard sphere - its great circle. So, the dynamics 
of the particle can be approximately described by how the position of this great circle slowly changes in time. 

o. 

Here we show that the result holds for any small smooth deformation of the sphere and that the procedure of 
obtaining the Hamiltonian of the reduced system consists in applying the Funk-Minkowski- Radon transformation (see 
of functions on the 2D-sphere to the function defining the difference of the surface from the sphere. 
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II. HAMILTONIAN REDUCTION 
A. Dynamics equations 

We consider the motion of a particle on a surface defined by the equation 

<p(x) = 

where 

3 

ip(x) = ^2x 2 t - 1 +eip(x), e«0. 
So, the surface is obtained by a small perturbation ip{x) of the standard sphere. 



FIG. 1: Two coils I, II of a trajectory on the surface (p(x) = 0; vectors Ni and Nu are the normals to the planes of great circles 
approximating the coils 



The motion of a particle with unit mass on the surface can be described by the Lagrange equations of the first kind: 



The Lagrange multiplier A can be found explicitly 



x ■ ——r ■ x x + e x ■ — — - ■ x 



\ = dxP _ = __ dxi_ (2) 

1 1 x + e 



dx J \ dx 
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Here -— 7 is the matrix of second derivatives of ib(x): 
ox 1 



fd 2 ip\ _ d 2 j; 
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B. Asymptotic descirption in terms of the momentum 



The trajectories of the particle have the following shape: small segments turn around the surface close to a plane 
section of it - a great circle (fig. [T]) . After a period of time this great circle changes its position on the surface. So, we 
can approximately describe the motion of the particle by tracing the slow change of the position of this great circle, 
thus determining at every moment of time the approximate location of the current loop of the geodesic. 

The position of a great circle can be described by a normal vector to it. It is easy to see that the angular momentum 
of the particle 

L = x x x 

provides such a normal vector. Its time derivative reads 

't - ^ - < dp - - r- , _ di\) ,„. 

L = x x x = x x A — = Xx x (x + e— -) = Ae x x — -. (3) 
ox ox ox 

Taking into account here only the first order terms in e, we can substitute A with its zero order component obtained 
from ©: Ao = —x 2 = const (the constancy follows from conservation of the kinetic energy E = x 2 /2). In what 
follows we assume the units to be chosen so that x 2 — 1 and hence Ao = — 1. 
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C. Averaging the momentum equations 



We can take advantage of the fact that the system is close to an integrable one - geodesies on the standard sphere. 
This leads to the system having a fast variable - the phase of the motion along the current loop (close to a great 
circle), and slow variables defining the position of the loop. The slow variables (components of the momentum) are 
first integrals for the unperturbed system. So, we are in a position to apply the method of averaging, 0. To derive 
an asymptotic description of the motion of the particle in terms of the evolution of its angular momentum L, we 
average the above equation over a period of the basic solution (for the unperturbed sphere) with a given (constant) 
momentum L: 

Xjf(t) = cost e\(L) + s'mt e 2 (Z), 

where vectors ei (L) , e*2 (£) form an orthonormal basis together with e^(L) — L/L (the exact choice of ei(L), e^Z) 
does not affect the result of averaging). This is the motion along the great circle perpendicular to the momentum L. 

So, we can see that the above averaging is equivalent to the integration of a function on the sphere over great circles, 
i.e. the Funk transform (see |4||): 

(Fg(x))(L) = J g (cost ei(£) + sint e 2 (£)) dt, 

where ei(L) and c-iilS) are the above-mentioned unit vectors orthogonal to L and to each other (fig. [5]). 
The Funk transform has the following basic properties: 



1. Linearity: F(Xg + fih) = \F(g) + fiF(h). 
This follows from the linearity of integration. 

2. Evenness of the image functions: (F g)(—L) — (F g)(L) (so, the image of the transform is a function on the 
projective plane). 

This is due to the fact that L and — L are perpendicular to the same great circle. 

3. Zero image of odd functions: if g(—x) = —g(x) then F g = 0. 
Indeed, for such g(x) we have 

(F g(x))(L) = Jq g (cost ei(Z) + sin t e 2 (L)j dt + g (cost ei(Z) + sin t e 2 (L)j dt — 
— Jq g (cost e\(L) + smt C2{L)\ dt + J™ g ( cos(u + n) e\(L) + sin(u + n) e2(L)j du = 
= Jq g { cost ei(L) + s'mt e2(L) ) dt + Jq g [ — cos it ei(L) — sinu e2(L) ) du = 0. 



So, averaging of the equation for the momentum ([3]) gives 

- d 

where we have denoted the operator i = ix — . 

ox 



D. The key commutation relation for the Funk transform 



Now we are going to derive the following general identity: 

F o f = To F, 



(5) 
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FIG. 2: The Funk transform of a function g(x) on the unit sphere, taken at the point L of the sphere, is the integral of g over 
a great circle perpendicular to L. 



where I is the operator 



= xx 



For that end we employ the following property of the Funk transform. 



E. Commutation of the Funk transform with rotations 

The Funk transform commutes with the action of three-dimensional rotations. Indeed, for any rotation R G 50(3) 
acting on functions as 

Rip(x) = tpiR^x) 

we have 

{R(Ftp))(L) = {Fip){R- l L) = [ ip{x)ds= [ ip{R~ l Rx) ds = [ ip{R~ l y) ds = {F R^){L), (6) 

i.e., we have obtained that 

Ro F ~ F o R. (7) 



F. Proof of the key commutation relation 



Now we note that the components of the operator I are the infinitesimal rotations: 



\a=0, 
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where Ri(a) is the rotation about the axis with number i through the angle a. E.g., for i = 3 we have the rotation 
matrix 

cos a — sin a ' 
Rs(a) = | sin a cos a 
1, 



and hence 



— (R3(a)ip(x)) = — ^(R 3 {a) 1 x) — — ih(x\ cosa + X2 sin a, —x\ sin a + X2 cos a, £3) 
ao aa da 

= —(xi cos a + X2 sin a)— h (xa cos a — xi sin a)— — 

0x2 oxi 



At a = we obtain 



-(R 3 (aMx))\ a=0 = - Xl — +X2 — = -lxx-) = 



3 

So, in ([7]), substituting i? = i?j(a), taking the derivative with respect to a, and substituting a = 0, we obtain the 
commutation relation ([5]): 

ft oF)i> = -^- ((Ri(a) o F) ^) | Q=0 = ~ ((F o ifc(a)) V) |a=o = (-F ° I<) ^ 
da aa 



G. Hamiltonian form of the momentum equations 

This observation enables us to change the places of operators in Q : 

t=-±-e (toF1,)=-±;e (Zx°)(FM£). (8) 



dL 

But it is easy to see that L x -Q* acts on a given function G(L) as taking its Poisson bracket with L: 

3C 

Lx- ? = {G,L}, 
oL 

where the usual Poisson brackets for the components of the momentum are assumed: 

{Li,Lj} — S ejjkLk- 
k 

So, we obtain that the averaged equations for the momentum can be formulated in the following Hamiltonian form: 

L = {L,H(L)}, (9) 
where the Hamiltonian is the Funk transform of the deformation function: 

H{L) = ±-F{ei,). 
III. CONCLUSION 

The reduced Hamiltonian system has L 2 as its first integral (a Casimir function), so the phase space is the 2D- 
sphere. Being a Hamiltonian system with one degree of freedom, it is completely integrable. The topological structure 
of the system's dynamics can be studied by means of phase portraits on the sphere L 2 = const. The trajectories are 
the contour lines of the Hamiltonian. 
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It is important to note that, in the generic case when the Hamiltonian is a Morse function, the topological char- 
acterization of the phase portrait can be performed using the topological classification of Morse functions on the 
2D-sphere developed by V.I. Arnold in Q. 
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